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Abstract 

We study the max-algebraic analogue of equations involving ^-matrices 
and M-matrices, with an outlook to a more general algebraic setting. We 
show that these equations can be solved using the Frobenius trace down 
method in a way similar to that in non-negative linear algebra 113], 1211 133] , 
characterizing the solvability in terms of supports and access relations. We 
give a description of the solution set as combination of the least solution 
and the eigenspace of the matrix, and provide a general algebraic setting 
in which this result holds. 

Keywords: Max-algebra; nonnegative linear algebra; idempotent semir- 
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1 Introduction 



A Z-matrix is a square matrix of the form XI — A where A is real and A is 
an (elementwise) nonnegative matrix. It is called an M-matrix if A > p (A) , 
where p (A) is the Perron root (spectral radius) of A and it is nonsingular if and 
only if A > p (A) . Since their introduction by Ostrowski [30] M-matrices have 
been studied in many papers and they have found many applications. The term 
Z-matrix was introduced by Fiedler-Ptak [TJ] . 

Results on the existence, uniqueness and nonnegativity of a solution x of 
the equation (XI — A) x = b for a given nonnegative vector b appear in many 
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places (e.g. Berman-Plemmons [1] in the case of a nonsingular M-matrix or 
an irreducible singular M-matrix). Using the Frobenius normal form of A and 
access relation defined by the graph of the matrix, Carlson [5] studied the exis- 
tence and uniqueness of nonnegative solutions x of this equation in the case of a 
reducible singular M-matrix, and his results were generalized to all Z-matrices 
in Hershkowitz-Schneider |21j . 

The purpose of the current paper is to prove corresponding results in the 
max times algebra of nonnegative matrices, unifying and comparing them with 
the results in the classical nonnegative linear algebra. We also notice that the 
basic proof techniques are much more general. In particular, we exploit a gen- 
eralization of the Frobenius trace down method [T3J [33J . This generalization is 
reminiscent of the universal algorithms developed by Litvinov, Maslov, Maslova 
and Sobolevskh [27J [28j [29] , based on the earlier works on regular algebra ap- 
plied to path-finding problems by Backhouse, Carre and Rote [UGH]. Following 
this line allows to include other examples of idempotent semirings, such as max- 
min algebra |16j and distributive lattices |35| . A more general theoretic setup is 
described in Section 0] It is very close to Cohen, Gaubert, Quadrat, Singer [DJ 
and Litvinov, Maslov, Shpiz [26 . 

The main object of our study is Z-matrix equation 

Ax + b = Xx. (1) 

over semirings. In the classical nonnegative algebra and max-plus algebra, any 
A 7^ is invertible, which allows to reduce ([1]) to 

Ax + b = x. (2) 

In max-plus algebra, this equation is sometimes referred to as discrete Bellman 
equation, being related to the Bellman optimality principle and dynamic pro- 
gramming [TJ 1181 122] . In particular, it is very well-known that this equation 
has the least solution. However (to the authors' knowledge) a universal and 
complete description of solutions of Xx = Ax + b or even x — Ax + b, which 
would cover both classical nonnegative and max-plus algebra cases, is not found 
(surprisingly) in the key monographs on max-plus algebra and related semirings. 
Such a description is what we try to achieve in this paper, see Theorems 13.21 
and 13.81 In brief, the results in the case of max times linear algebra are similar 
to those in classical Z- matrix theory |21j , but they are not identical with them. 
Details are given in the main sections. The situation is analogous to that for 
the Perron-Frobenius equation Ax = Xx, as may be seen by comparing the re- 
sults in Butkovic [7], Sect. 4.5, for max algebra with those in Hershkowitz and 
Schneider [2D], Sect. 3, for classical nonnegative algebra. 

The rest of the paper consists of Prerequisites (Sect. 2), Theory of Z-matrix 
equations (Sect. 3) and Algebraic generalization (Sect. 4). Prerequisites are 
devoted to the general material, mostly about max-plus algebra: Kleene star, 
Frobenius normal forms and spectral theory in the general (reducible) case. 
Theory of Z-matrix equations over max-plus algebra stands on two main re- 
sults. Theorem 13.21 describes the solutions of ([2]) as combinations of the least 
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solution A* b and the eigenvector space. We emphasize the algebraic generality 
of the argument. Theorem 13.81 exploits the Frobenius trace down method. This 
method serves both for theoretic purposes (to provide a necessary and suffi- 
cient condition for the existence of solutions, and to characterize the support 
of the least solution) and as an algorithm for calculating it. As an outcome, 
we get both combinatorial and geometric description of the solutions. The re- 
sults in max-times algebra are compared with the case of nonnegative matrix 
algebra [21] ■ The paper ends with Section U] devoted to an abstract algebraic 
setting for which Theorem [321 holds, in the framework of semirings, distributive 
lattices and lattice-ordered groups [5J [T7] . 

We use the conventional arithmetic notation a+b and ab for the operations in 
semirings, to emphasize the common features of the problem over classical non- 
negative algebra and in the idempotent semirings, viewing max-times algebra 
(isomorphic to max-plus algebra) as our leading example. 

We note that a complete description of solutions of x = Ax + b was also 
achieved by Krivulin [231 |2~31 |2"5] , for the case of max algebra and related semir- 
ings with idempotent addition (a © a = a) . His proof of Theorem 13.21 recalled 
here in a remark following that theorem, is different from the one found by the 
authors. We show that Krivulin's proof also works both for max algebra and 
nonnegative linear algebra, and admits further algebraic generalizations. The 
case of reducible matrix A and general support of b has been also investigated, 
see [53] Theorem 2 or [55] Theorem 3.2, which can be seen as a corollary of The- 
orem 13.81 of the present paper with application of the max-plus spectral theory, 
see Theorem 12.21 

2 Prerequisites 

2.1 Kleene star and the optimal path problem 

The main motivation of this paper is to unite and compare the Z-equation 
theory in the classical nonnegative linear algebra, and the max-times linear 
algebra. Algebraically, these structures are semirings |17] (roughly speaking, 
" rings without subtraction" , see Section [4] for a rigorous definition) . Thus we 
are focused on 

Example 1: Max-times algebra. Nonnegative numbers, endowed with 
the usual multiplication x and the unusual addition a + b := max(a,6). 

Example 2: Usual nonnegative algebra. Nonnegative numbers, en- 
dowed with usual arithmetics +, x. 

Some results in this paper will have a higher level of generality, which we 
indicate by formulating them in terms of a " semiring 5" . Namely, this symbol 
"5" applies to a more general algebraic setting provided in Section [4[ covering 
the max-times algebra and the nonnegative algebra. Before reading the last 
section, it may be assumed by the reader that S means just "max-times or 
usual nonnegative" . 

The matrix algebra over a semiring S is defined in the usual way, by extend- 
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ing the arithmetical operations to matrices and vectors, so that (A + B)ij = 
a ij + and (AB)ik = J2j a ijbjk for matrices A,B of appropriate sizes. The 
unit matrix (with l's on the diagonal and O's off the diagonal) plays the usual 
role. 

Denote N = {1, n} . For x £ S n , x > means Xi > for every i. Similarly 
A > for A e S nxn . We also denote: 

A* =I + A + A 2 + ... = sup(7 + A + ...+ A k ). (S) 

k>0 y ' 

In we have exploited the nondecreasing property of addition. A* is also 
called the Kleene star, and it is related to the optimal path problem in the 
following way. 

The digraph associated with A = (aij) £ S nxn is Da = (N,E), where 
E = {(i,j) ; > 0}. The weight of a path on Da is defined as the product of 
the weights of the arcs, i.e., the corresponding matrix entries. It is easy to check 
(using the distributivity law) that (A k )ij is the sum of the weights of all paths 
of length k connecting i to j. Further, an entry (A*)ij collects in a common 
summation (possibly divergent and formal) all weights of the paths connecting 
i to j, when i j. 

Note that A* = (I — A)~ x in the case of the classical arithmetics, and A* 
solves the optimal path problem in the case of the max-times algebra (because 
the summation is maximum). 

Thus the Kleene star can be described in terms of paths or access relations 
in Da- For i =/= j, we say that i accesses j, denoted i — > j, if there is a path 
of nonzero weight connecting i to j, equivalcntly, (A*)ij ^ 0. We also postulate 
that i — > i. The notion of access is extended to subsets of N, namely, / — > J 
if i — > j for some i € I and j G J. 

Both in max-plus algebra and in nonnegative algebra the Kleene star se- 
ries may diverge to infinity (in other words, be unbounded). In both cases the 
convergence is strongly related to the largest eigenvalue of A (w.r.t. the eigen- 
problem Ax = \x), which we denote by p(A). This is also called the Perron 
root of A. A necessary and sufficient condition for the convergence is p(A) < 1 
in the case of the ordinary algebra, and p(A) < 1 in the case of the max-times 
algebra. In the max-times algebra (but not in the usual algebra) A* can be 
always truncated, meaning A* = I + A + . . . + A" -1 where n is the dimension 
of A, in the case of convergence. This is due to the finiteness of Da and the 
optimal path interpretation, see p] [7] for more details. 

In the case of max-times algebra, p(A) is equal to the maximum geometric 
cycle mean of A, namely 

p(A) = max { tyai 1 i 2 ai 2 i 3 ...a ik i 1 ; i%, i k E N,k = 1,2,...} . 

This quantity can be computed in 0(n 3 ) time by Karp's algorithm, see e.g. [H[7]. 

2.2 Probenius normal form 

A = (a^) € K" x ™ is called irreducible if n = 1 or for any i,j £ N there 
are i\ = i,i 2 ,---,ik = j, such that 

&Z1Z2 i 3 "- a i k -ii k > 0', A is called reducible 
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otherwise. In other words, a matrix is called irreducible if the associated graph 
is strongly connected. Note that if n > 1 and A is irreducible then A ^ 0. Hence 
the assumption "A irreducible, A ^ 0" merely means that A is irreducible but 
not the lxl zero matrix. (It is possible to extend these notions to general 
semirings with no zero divisors, but we will not require this in the paper.) 

In order to treat the reducible case for max-times algebra and the (classical) 
nonnegative linear algebra, we recall some standard notation and the Frobenius 
normal form (considering it for general semirings will be of no use here). If 

1 < h < i 2 < ... < ik < n,K= {i!,...,i k } C N 

then Akk denotes the principal submatrix 



of the matrix A = (ay) and xk denotes the subvector [xi x , Xi k ) T of the 
vector x = (ii, x n ) T . 

If D = (N,E) is a digraph and K C N then D(K) denotes the induced 
subgraph of D, that is 

D(K) — (K, E (1 (K x K)). 

Observe that p (A) = if and only if Da is acyclic. 

Every matrix A — (ay) G R!f: x ™ can be transformed by simultaneous per- 
mutations of the rows and columns in linear time to a Frobenius normal form 

ED 

( A n ... \ 

A 2 i A 22 ... ^ 

y A r i A r2 ... A rr f 

where An, A rr are irreducible square submatrices of A, corresponding to the 
partition Ni U . . . U N r — N (that is, Ay is a shortcut for A^at )• The sets 
Nx,...,N r will be called classes (of A). It follows that each of the induced 
subgraphs DA(Ni) (i = 1, r) is strongly connected and an arc from Ni to Nj 
in Da exists only if i > j. 

If A is in the Frobenius normal form Q then the reduced graph, denoted 
R(A), is the digraph whose nodes are the classes Ni,...,N r and the set of arcs 
is 

{{N U N ); (3k e NA{31 E N 3 )a M > 0}). 

In addition we postulate that each class has a self-loop (useful if FNF contains 
trivial classes consisting of one diagonal zero entry). In the max-times algebra 
and the nonnegative matrix algebra, the nodes of R(A) are marked by the 
corresponding greatest eigenvalues (Perron roots) pi := p(Au). 

Simultaneous permutations of the rows and columns of A are equivalent 
to calculating P~ 1 AP, where P is a generalized permutation matrix. Such 
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transformations do not change the eigenvalues, and the eigenvectors before and 
after such a transformation only differ by the order of their components. So 
when solving the eigenproblem, we may assume without loss of generality that 
A is in a Frobenius normal form, say Q. 

2.3 Eigenvalues and eigenvectors in max-times algebra 

It is intuitively clear that all eigenvalues of A are among the unique eigenvalues 
of diagonal blocks. However, not all of these eigenvalues are also eigenvalues 
of A. The following key result, describing the set A(A) of all eigenvalues of A 
in max-times algebra (that is, set of all A such that Ax — Xx has a nonzero 
solution x in max-times algebra) appeared for the first time independently in 
Gaubert's thesis [13] and Bapat et al. [3], see also Butkovic [7]. 

Theorem 2.1 ((cf. [3l W \ fT4] ^ Let ^ be a Frobenius normal form of a ma- 
trix A e M+ Xn . Then 

A (A) = {pj;pj > Pl for all JV< — -> Nj}. (5) 

The same result holds in the nonnegative linear algebra, with the non-strict 
inequality replaced by the strict one. 

If a diagonal block Ajj has pj G A, it still may not satisfy the condition in 
Theorem 12.11 and may therefore not provide any eigenvectors. So it is necessary 
to identify classes j that satisfy this condition and call them spectral. Thus 
Pj € A(A) if Nj is spectral, but not necessarily the other way round. We can 
immediately deduce that all initial blocks are spectral, like in the nonnegative 
linear algebra. Also, it follows that the number of eigenvalues does not exceed 
n and obviously, p(A) = max^p^, in accordance with p(A) being the greatest 
eigenvalue. 

We are now going to describe, for A € A, the eigencone V(A, A) of all vectors 
x such that Ax — Xx. Denote by J\ the union of all classes Ni which have access 
to the spectral classes corresponding to this eigenvalue. By (0, p.; < A for all 
such classes. Now we define the critical graph Ca(X) = (N C ,E C ) comprising all 
nodes and edges on critical cycles of the submatrix Aj x j x , i.e., such cycles where 
A is attained. This graph consists of several strongly connected components, 
and let T(A, A) denote a set of indices containing precisely one index from each 
component of Ca(X). In the following, A'(J\) will denote the n x n matrix, 
which has Aj x j x /X as submatrix, and zero entries everywhere else. 

Theorem 2.2 ((cf. HO [14])) Let A e K" xn and X e A(A). Then 

a) For any eigenvector v € V{A, X) there exist otj e M + such that v is the 
max-times linear combination 

v= ]T ajiA'iJxWj. (6) 

3 £T(A,\) 
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b) For any two indices j and k in the same component of CU(A), columns 
(A'(J\))*j and (A'(J\))* k are proportional. 

c) Vectors (A'(J\))* for j g T(A, A) form a basis ofV(A,\), that is, they 
generate V(A, A) in the sense of a) and none of them can be expressed as 
a max-times linear combination of the others. 

Remark 2.3 An analogous description of V(A, A) in nonnegative matrix alge- 
bra is called Frobenius- Victory theorem (T3J[T!5], see [33J Theorem 3.7. Namely, 
to each spectral node of R(A) with eigenvalue A, there corresponds a unique 
eigenvector, with support equal to the union of classes having access to the 
spectral node. These eigenvectors are the extreme rays of the cone, i.e. they 
form a "basis" in analogy with Theorem 12.21 

Moreover, these extreme rays are linearly independent as it may be deduced 
from their supports. In parallel, it can be shown that the generators of Theo- 
rem !2.2l are strongly regular, see [5] for definition (i.e., independent in a stronger 
sense) . 

However, extremals in the nonnegative case do not come from A* = (/—A) -1 
and, to the authors' knowledge, no explicit algebraic expression for these vectors 
is known. 

3 Theory of Z-matrix equations 
3.1 General results 

In the following, we describe the solution set of x = Ax + b as combinations of 
the least solution A*b and the eigenvectors of A. The results of this subsection 
hold for the max-times algebra, nonnegative linear algebra, and an abstract 
algebraic setup (reassuring that A*b satisfies x = Ax + b, and v = inf A k x, for 

k 

x such that Ax < x, satisfies Av = v), which will be provided in Section |4l 
We start with a well-known fact, that 

A*b := sup(6 + Ab + A 2 b + ... + A^b) (7) 

k 

is the least solution to x = Ax + b. We will formulate it in the form of an 
equivalence. Note that the supremum ((7|) may exist even if A* does not exist. (In 
this sense, A*b is rather a symbol than a result of matrix- vector multiplication. 
On the other hand, one can complete a semiring with the greatest element " +oo" 
and regard A* b as a matrix- vector product.) 

Theorem 3.1 ((Well-known, cf. [HUH])) Let A € S nxn ,b € S n . The fol- 
lowing are equivalent: 

(i) x — Ax + b has a solution, 

(ii) x = Ax + b has a least solution. 
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(Hi) A*b converges. 

If any of the equivalent statements holds, A*b is the least solution of x = Ax + b. 

Proof. (i)=> (iii) Let a; be a solution to x = Ax + b. Then 

x = Ax + b 

= A (Ax + b) + b 

= A(A(Ax + b) + b) + b= .... 

Therefore for any k > 1 we have 

x = A k x+ (A 1 "- 1 + A k ~ 2 + ... + I)b. (8) 

This shows that the expressions in ([7]) are bounded from above by x, hence the 
supremum exists. 

(iii)=> (ii) We verify that 

AA*b + b = A sup(6 + Ab + . . . + i4* _1 6) + b 

k v ( 9 ) 

= sup(6 + Ab + . . . + A k b) = A*b, 

k 

treating sup as a limit and using the continuity of the matrix- vector multiplica- 
tion. 

(From the algebraic point of view, we used the distributivity of (max- algebraic, 
nonnegative) matrix multiplication with respect to sup's of ascending chains, 
and the distributivity of + with respect to such sup's. Further details on this 
will be given in Sectional) 

(ii)=*> (i) Trivial. ■ 

We proceed with characterizing the whole set of solutions. (See also Re- 
mark [33] for an alternative short proof of the first part.) 

Theorem 3.2 Let A E S nxn ,b £ S" be such that x = Ax + b has a solution. 
Then 

(a) The set of all solutions to x — Ax + b is {v + A*b; Av = v}; 

(b) If for any x such that Ax < x we have inf A k x — 0, then A*b is the unique 
solution to x — Ax + b. 

Proof, (a) Firstly we need to verify that any vector of the form v + A*b where 
v satisfies Av — v, solves @. Indeed, 

A{v + A*b) + b = Av + (AA*b + b) = v + A*b, 

where we used that Av = v and that A*b is a solution of @, see Theorem 13. II 
It remains to deduce that each solution of is as defined above. 
Let a; be a solution to ([2J, and denote ?/ fe ' := A k x and 

z (k) ,_ (a**- 1 + A k - 2 + ... + I) b. 
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We have seen in © that 

x = y {k) + z (k \ forallfc>l. (10) 

Since Ax < x it follows that the sequence y^ k ' is nonincreasing. The sequence 
of z^ is nondecreasing. 

Both in max-times and in the nonnegative case, we conclude that v = 
lim/^oo z/ fe ) exists and (by the continuity of A as operator) we have Av = v. 
We also obtain that A*b — limfc_ s . 00 z^ k \ and finally 

x = lim y {k) + lim z {k) =v + A*b, (11) 

where i; satisfies Au = v. The theorem is proved, both for max-times algebra 
and nonnegative linear algebra. 

In a more general semiring context (see Section^), it remains to show that 
x, := inffe y( k > + sup fe z^ is the same as y^ + z^ for all k. After showing this 
we are done, since supz( fc ) = A*b, and also 

k 

A inf y (k) = Afinf A k x) = inf A k x = inf y {k) , (12) 

fe>0 fc>0 k>l k>0 

so that we can set v := inffc y^ k \ it satisfies Av — v. (From the algebraic point 
of view, we have used the distributivity of matrix multiplication with respect to 
inf's of descending chains. Further details will be given in Section |4j) 
Using the distributivity of + with respect to inf we obtain 

x = inf(y (fe) +supz W ) > x, (13) 
k I 

since this is true of any term in the brackets. Using the distributivity with 
respect to sup we obtain 

x = sup(inf y {k) + z {l) ) < x, (14) 
; k 

for analogous reason. Combining (jTBf and (|T4")) we obtain 

x = x = inf j/ fe) + sup z (fe) = v + A*b, 

k k 

which yields a general proof of part (a). 

For part (b), recall that ?/ fe ) :— A k x, and that x satisfies Ax < x. ■ 
These results also apply to equations Xx — Ax + b when A is invcrtiblc: it 

suffices to divide this equation by A. 

Remark 3.3 The solution set of x = Ax + b is convex over S, since it contains 
with any two points x, y all convex combinations Xx+fj,y, A+/i = 1. Further, both 
in max-times semiring and in the nonnegative algebra, A*b is the only extreme 
point: it cannot be a convex combination of two solutions different from it. 
The eigenvectors of A are recessive rays,i.e., any multiple of such vectors can 
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be added to any solution, and the result will be a solution again. Moreover, 
only eigenvectors have this property. Indeed, assume that z is a solution, z + fiv 
where fj, =^ satisfies 

z + fiv — A(z + jjLv) + b, 

but Av 7^ v. In the usual algebra this is impossible. In max-times, assume 
that (Av)i ^ Vi for some i, then one of these is nonzero. As zi = (Az)i + bi 
is finite, taking large enough /i will make fiVi or (j,(Av)i the only maximum 
on both l.h.s. and r.h.s., in which case z + /ad will not be a solution. Thus, 
in both theories the eigencone of A with eigenvalue 1 is the recessive cone of 
the solution set of x = Ax + b. In the max-times case it is generated by the 
fundamental eigenvectors as in Theorem 12.21 Thus we have an example of the 
tropical theorem of Minkowski, representing closed max-times convex sets in 
terms of extremal points and recessive rays, as proved by Gaubert and Katz [15] . 

Remark 3.4 In this remark we recall the proof of Theorem [32] part a) given by 
Krivulin, see [23] Lemma 7 or [25] Lemma 3.5. We slightly change the original 
proof to make it work also for nonnegative linear algebra. Let s be a solution 
of x = Ax + b and define w as the least vector w satisfying x = u + w where 
u := A*b. It can be defined explicitly by 

{•Eii if <Ei ^ 'Wj: I Xi Uj, if Xl ^> V>i^ . . 

„ ., or w t = ^ , (15) 

0, if x l = m, [0, it Xi = Ui, 

in the case of max algebra and nonnegative linear algebra respectively. Now 
notice that if x — u + w then x — u + Aw. Indeed 

x = A(u + w) + b = (Au + b) + Aw = u + Aw. 

Hence w < Aw. Indeed, both w := w and w := Aw satisfy x = u + w but w 
is the least such vector. Defining v := sup n>1 A n w we obtain x = u + v and 
Av = v. The algebraic generality of this argument is also quite high, it will be 
discussed in the last section of the paper. 



3.2 Spectral condition, Probenius trace-down method 

We consider equation Xx = Ax + b in max-times algebra and nonnegative 
linear algebra, starting with the case when A is irreducible. Theorem 13.81 
below can be also viewed in nonnegative linear algebra, but only after some 
modification which will be described. Denote A\ := A/X and b\ := b/X. 

The following is a max-times version of the Collatz-Wielandt identity in 
the Perron-Frobenius theory. 

Lemma 3.5 ((Well-known, cf. [3 HI])) Let A G R™ x ", A ^ 0. Then Ax < 
Xx has a solution x > if and only if A > p (A) , A > 0. 
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Proof. Let x > be a solution, then Ax ^ and so A > 0. II p (A) = there 
is nothing to prove, so we may suppose p (A) > 0. Let a — ik, ik+i = h) 

be any cycle with nonzero weight. Then 

^122*^12 — Ax^ 
^1213*^3 — AXj 2 

^ikii^ii — XXi k . 

After multiplying out and simplification we get A > tya,i 1 i 2 ai 2 i 3 ...a,i k i 1 and so 
\>(>(A). 

Suppose now A > p (A) , A > 0. Then p (A\) < 1 and so A* x = I+A\ + ...+A k x 
for every k > n — 1, yielding AaA^ < A^. Let it be any positive vector in R™ . 
Take x = A x u, then x > because ^4^u > u and 

A A a; = A x A\u < A* x u = x. 



Lemma 3.6 ((Well-known, cf. [H[30])) If A e M" x " is irreducible, b € 
WL, b and A > then the following are equivalent: 

(i) Xx — Ax + b has a solution. 

(ii) A* x converges. 

(Hi) A > p(A) (max-times algebra), A > p(A) (nonnegative linear algebra). 
All solutions of \x = Ax + b (if any) are positive. 

Proof. In the case of nonnegative matrix algebra, this Lemma follows from 
the results of Ostrowski's famous paper [30] (where matrices appear as deter- 
minants), see also [32], Lemma 5. For the equivalence between (ii) and (iii) in 
max-times algebra, consult e.g. [J [7] US] , (Both in max-times algebra and in 
the nonnegative linear algebra, such equivalence holds also for reducible matri- 
ces.) For the reader's convenience we show the equivalence between (i) and (iii) 
in max-times algebra. 

(iii)=>(i): If A > p (A) then 1 > p(A\), hence A* x and A* x b\ converge. In this 
case, A* x b\ is the least solution by Theorem 13. II 

(i)=>(iii): If Ax + b = Xx then A > and x ^ since b ^ 0. We need to show 
that x > 0, to apply Lemma l3~5l 

If n = 1 then the result holds. Suppose now n > 1, thus p (A) > 0. Let 
B = {piA))^ 1 A and p = (p(A)) -1 A. Then B has p(B) = 1, it is irreducible 
and Bx < px. Therefore B* > 0, thus B*x > 0. But B*x < px and hence 
x > 0. By Lemma T3. 5 1 we obtain that A > p(A). ■ 

Remark 3.7 Note that also for general (reducible) A, if b > then for any 
solution x of Ax = Ax + b we have x > 0, and hence A > p (A) by Lemma 13.61 
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However, this condition is not necessary for the existence of a solution to Ax + 
b = Xx, when b has at least one zero component, see Theorem 13.81 below. If 
A < p(A) then some entries of A* x are +co and it is not obvious from Theorem 
13.21 whether a finite solution exists since the product A* x b\ may in general (if A 
is too low) contain +00. However if 0. (+00) is defined as and the condition 
of Theorem 13.81 (iv) holds, then the +00 entries of A* x in A\b\ will always be 
matched with zero components of b\, and consequently A x b\ will be a finite 
non- negative vector. 

Now we consider the general (reducible) case. The next result appears as 
the main result of this paper, describing the solution sets to Z-equations in 
max-times algebra and nonnegative linear algebra. 

Theorem 3.8 Let A e K" x " be in FNF with classes Nj,j = l,...,s. Let 
b £ R™,A > 0,. Denote J = {j',Nj — !>supp(6)} and ~p = max je j pj (for 
the case when b — and J = assume that max0 = 0). The following are 
equivalent: 

(i) System Xx = Ax + b has a solution. 

(ii) System Xx = Ax + b has the least solution. 
(Hi) x^ = A* x b\ converges. 

(iv) If j € J then (Ajj) x converges. 

(v) ~p < X (max-times), or~p < A (nonnegative linear algebra) 
If any of the equivalent statements hold, then 

a) x° is the least solution of Xx = Ax + b. For this solution, x^, 7^ when 
i £ J and x^. — when i ^ J . The solution x° is unique if and only if X 
is not an eigenvalue of A. 

b) Any solution x of (Qp can be expressed as x = x° + v where v satisfies 
Av = Xv. 

Proof. We first treat the trivial case, when 6 = 0. In this case x° = is a 
solution, J = 0,p = O<A and thus all the equivalent statements (i)-(iv) are 
true; (a) and (b) hold trivially with x° = 0. 

We now suppose b 7^ 0. Consequently, A > 0, and assume w.l.o.g. that 
A = 1. The equivalence of (i)-(iii) was manifested in Theorem 13.11 and part 
b) was proved in Theorem 13.21 The equivalence of (iv) and (v) follows from 
Lemma l3.6l It remains to show the equivalence of (i) and (iv), that the minimal 
solution has a prescribed support, and the spectral criterion for uniqueness. 

We show that (i) implies (iv) . For simplicity we use the same symbol " J" 
for the set of indices in the classes of J. Denote / := {1, . . . , n}\ J. We have 
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and hence xj is a solution of Ajjxi = xj, and xj is a solution of xj — Ajjxj + 
AjiXi + bj. Further, denote bj := Ajjxj + bj and let J consist of the classes 
JVi,...,JV t (t< s). Then 



A(J) 



( A n 
A 2i v4 22 



A, 



\ 




A u j 



We now proceed by an inductive argument, showing that (Ajj)* converges 
for all j = 1, . . . , t, and that all components in xn x , . . . , xjy t are positive. This 
argument is a max-algebraic version of the Frobenius trace-down method. 

As the base of induction, we have xn x — AuXni + b^. In this case, the 
class N\ is final, so b]y 1 and hence bjvi should have some positive components. 
Using Lemma T3. 6 1 we conclude that (An)* converges and xn 1 is positive. 

Induction step. Suppose that for some I, all components of (xn^ ■ ■ ■ ,xn,) 
solving 



XN t = Anx Nl + b Nl 

xn 2 = A 2 ix Nl + A 2 2Xn 2 + b N . 2 

xn, = A a x Nl + . . . + AiiXn, + b Nl , 
are positive. We show that the same holds if we add the next equation 



(17) 



x N , 



Ai+i^XNt + • • . + A l+u+1 x Nl+1 + b 



n,. 



(18) 



and that (Ai + i } i + i)* converges. We have two cases: either bN l+1 has nonzero 
components so that N1+1 intersects with supp(fo), or if not, iV;+i should access 
a class which intersects with supp(6). In this case, one of the submatrices 
Ai+i^i, . . . , Ai+i t i is not identically 0. As all components of xn x , ■ ■ ■ ,xn, are 
positive, this shows that the sum on the r.h.s. of (1181) excluding Ai + ij + iXn 1+1 
has some positive components in any case. Using Lemma l3.6[ we conclude that 
(Ai + i t i + i)* converges and xn 1+1 is positive. 

Now we show that (iv) implies (i), and moreover, that there is a solution 
with prescribed support structure. To do so, we let xj = in (|16| . Then 
it is reduced to xj = Ajjxj + bj, and we have to show the existence of a 
positive solution xj. The proof of this follows the lines of the Frobenius trace- 
down method described above, making the inductive assumption that (|17p has 
a positive solution (xn x , ■ ■ ■ , xn, ) and using Lemma 13.61 to show that (|18|) can 
be solved with a positive xn 1+1 ■ Strictly speaking, in this case we have b instead 
of b in (|17|) and (|18|) . but this does not make any change in the argument. 

Let the conditions (i)-(v) be satisfied. Since letting xi = in P^)l produces a 
solution (see above), the support of the least solution is contained in J. However, 
the support of any solution should contain J by the argument in the proof of 

(i)=>(iv). 
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Evidently, solution x° is unique if A is not an eigenvalue of A. To show the 
converse (in max-times algebra), note that for any nonzero v there is a large 
enough a such that some component avi is greater than x®, hence x° + av ^ x° . 
(Note that the converse would be evident in the usual nonnegative algebra.) 

The proof is complete. ■ 

Remark 3.9 ((cf [2, 28j )) The Frobenius trace-down method of Theorem l3.8l 

can be also viewed as a generalized block-triangular elimination algorithm for 
obtaining the least solution A*b (assumed w.l.o.g. that A = 1). Namely, if 
(xn 1 . . . Xn, ) is the least solution of ([XT]) , then computing 

xn 1+1 ~ (Ai +lil+ i)*(A l+1:1 x Nl + . . . + A l+ ijx Nl + b Nl+1 ) (19) 

yields the least solution (xn x , . . . , xn { , %N l+1 ) of the enlarged system (|17l) & 
(fT5|) with b instead of b. Indeed, if we suppose that there is another solution 
(x' Ni , . . . ,x' Nii ), then x' N . > xjv 4 , and it follows from (fT9)) that x' Ni+i > xjy, +1 . 
As an algorithm for finding the least solution of x = Ax + 6, it is valid even for 
more general semirings than the setting of Section [4j provided that a solution 
to x = Ax + b exists. 

Remark 3.10 The description of the support of x° as in Theorem 13.81 a) can 
be obtained directly from the path interpretation of A*, using that x° = A*b 
(when this is finite). Indeed, write b = J2kesupp(b) Pu^k where is the fcth unit 
vector and j3k are all positive. Hence x° = A*b = X)fcesupp(b) PkA* k . It can 
be now deduced from the path interpretation of A*, that x® > whenever I 
accesses k from supp(6). This argument also shows that the description of the 
support of x° = A*b is valid over any semiring with no zero divisors. With zero 
divisors, the access condition for x® > may be no longer sufficient. 

We have represented any solution x of x = Ax + b in the form x° + v, where 
x o _ j^tjj an( j y _ j n £ fe A k x . Below we give an explicit formula for v in the case 
of max-times algebra, due to Dhingra and Gaubert |10) . 

Let C be the set of critical nodes (i.e., nodes of the critical graph) correspond- 
ing to the eigenvalue 1 (see Subsection l2.3[) . For any critical cycle («i, . . . , ik) we 
either obtain x^ = for all? = 1, . . . , k, or both x^ ^ and ai l i l+1 Xi l+1 — x^ for 
all I (after multiplying all inequalities ai t i l+1 Xi l+1 < x^ and canceling Xi t . . . Xi k 
it turns out that any strict inequality causes ai t i 2 . . ■a,i k i 1 < 1). This implies 
{Ax)c = xc, for the critical subvectors of Ax and x. Applying A to x^ — A^x 
which also satisfies Ax^ < x^ we obtain that (A k x)c = xc for any k, and 
hence also vc = xq- 

It remains to determine the non-critical part of v. For this we expand the 
non-critical part of Av — v as vn = A^cvc + A^nVn- Forming Aatjv corre- 
sponds to eliminating all spectral nodes of eigenvalue 1 from the reduced graph 
R{A). The non-spectral nodes with eigenvalue 1 will remain non-spectral, hence 
Ann does not have eigenvalue 1, and (Ann)* Anc^c is the only solution. Com- 
bining with the previous argument we conclude that 

v c = x c , v N = (Ann)* Ancxc- (20) 
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3.3 Max-times algebra and nonnegative linear algebra 

We can make further comparison with the survey of Schneider [33J Sect. 4, and 
with the work of Hcrshkowitz and Schneider |21] describing solutions of Z-matrix 
equations in nonnegative linear algebra. It can be seen that: 



1. In the nonnegative linear algebra, Theorem 13.81 extends the statements 
of [33J Theorem 4.3 and [33J Theorem 4.12 (due to Victory [19 ). In 
particular, it gives an explicit formula for the least solution. 

2. Frobenius trace-down method is also used in the proof of [21] Proposi- 
tion 3.6, and condition (v) of Theorem 13.81 (with the strict inequality) is 
equivalent to 21 condition (3.12). 

3. As observed in [21] Theorem 3.16, in the case of nonnegative algebra, x° 
is the only vector with support J, because supp(6) cannot be accessed by 
the spectral classes of R(A) in the case of solvability. However, this is 
not the case in the max-times algebra, making it possible that all spectral 
classes of R(A) access supp(fe). This condition is necessary and sufficient 
for all solutions of Xx = Ax + b to have the same support as A* b. 

4. It can be observed that geometric and combinatorial descriptions of the 
solution set in the usual nonnegative algebra, as provided by 2l[ Theo- 
rem 3.20 and Corollary 3.21, can be deduced from Theorem 13.81 with an 
application of Frobenius- Victory theorem, see remark after Theorem 12.21 
Max-times analogues of these results of [21] can be also easily formulated. 

We next give an example illustrating similarity and difference between the 
two theories. Let A be the matrix 

(l 0\ 

1 

1 1 
1 
1 
1 1 

\0 1 2/ 

We note that this matrix is essentially the same as in [3T] Example 3.22, that is, 
we have replaced / — A by A and its (reduced) graph R(A), given below, differs 
from the one in that Example 3.22 only by the markings of the nodes. 

Let b € M7 + . It follows from Theorem 13.81 (iv) that there exists a solution x 
to the max times equation Ax + b = x if and only if supp(&) C {1, 3, 4, 6}. In the 
usual nonnegative algebra, the condition is more restrictive: supp(6) C {4,6}. 

We choose 

b=(0 1 0) T 

as in [21] Example 3.22. Then supp(fe) = {4} and the minimal solution x° of 
Ax + b = x has support {4, 6} and equals 



x 



-° = (0 1 1 o) T . 
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1 1 

10 



Figure 1: The marked (reduced) graph of matrix A. Circles correspond to the 
nodes with the greatest Perron root 1, the darker ones being spectral. Each 
node is marked by its Perron root (inside) and by its number (outside) . 

in both theories. 

In max-times algebra, {1} and {3} are spectral nodes, and the eigenvector 
cone for the eigenvalue 1 is generated by 

v 1 = (1 1 1 0) T 

and 

v 2 = (0 1 1 o) T , 

see Theorems 12.11 and 12.21 In the usual nonnegative algebra, {3} is the only 
spectral node and any eigenvector is a multiple of v 2 . 

In max-times algebra, the maximal support of a solution is {1, 3, 4, 6}. For 
example take 

y 1 = (2 3 1 3 0) T , 

the max-times "sum" of x°, 2v x and 3v 2 . In the usual nonnegative algebra, the 
maximal support is {3, 4, 6}, take 

y 2 = x° + v 2 = (0 1 1 2 0) T , 

as in [21]. Note that neither y 1 is a solution in the usual sense, nor y 2 is a 
solution in the max-times sense. 

Observe that for given A, b, if the usual Perron roots of all blocks in FNF 
are the same as the max-times roots (as in the example above), the existence of 
a solution in nonnegative algebra implies the existence in max algebra (but not 
conversely) . Examples of vectors for which a solution exists with A as above in 
max algebra but not in nonnegative algebra are given by v 1 and v 2 . 

In the case of existence, minimal solutions have the same support in both 
theories. 
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4 Algebraic generalizations 



Here we describe general setting in which Theorems l3.ll and l3.2l of Subsection l3.ll 
hold. 

Recall that a set S is called a semiring if it is equipped with operations of 
addition + and multiplication • satisfying the following laws: 

(a) addition is commutative a + b = b + a Va, b G S, 

(b) multiplication distributes over addition a(b + c) = ab + ac Va, 6, c G S, 

(c) both addition and multiplication are associative: a+ (b + c) = (a + b) + c, 
a(6c) = (a6)c Va, o, c G 5, 

(d) there are elements and 1 such that a + = a, al = la = a, and 
aO = Oa = for all a G 5. 

The max-times algebra and the usual nonnegative algebra are semirings (also 
note that any ring and any field is a semiring), see also other examples below. 

We consider a semiring S endowed with a partial order <, i.e., binary relation 
< such that 1) a < b, b < c imply a < c, 2) a < 6, b < a imply a — 6, 3) a < a. 
In the case of idempotent addition (a + a = a), one defines a canonical order by 
a<&<=^a + & = 6. 

To model both max-times algebra and the usual nonnegative algebra, we 
may assume that the following axioms are satisfied. 

Al. any countable ascending chain (i.e., linearly ordered subset) in S bounded 
from above has supremum, and any countable descending chain has infi- 
mum; 

A2. addition is nondecreasing: a + b > a and a + b > b; 

A3, both operations distribute over any supremum or infimum of any such 
chain in 5, i.e. 

a + sup6 M = sup(a + & M ), a' • supo M = sup(a' ■ 6 M ) 

MM M M ^21) 

c + inf gL = inf (c + dn), c ■ inf gL = inf(c' • cL) 
mm' mm' 

for any countable bounded ascending chain {6 M } C S, countable descend- 
ing chain {d^} C 5, elements a, c, a', c' G S. 

Axiom A2 implies that the semiring is nonnegative: a > for all a, and 
antinegative: a+b — implies a = b = 0. Axiom A3 implies that both arithmetic 
operations are monotone. 

The operations of S extend to matrices and vectors in the usual way. More- 
over we can compute matrix powers A k for k > 0, where we assume A° = I, the 
identity matrix, where all diagonal entries equal to 1 and all off-diagonal entries 
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equal to 0. The extension of notions of associated digraph and access relations 
is also evident, provided that there are no zero divisors. 

Note that partial order in S is extended to S n and S mxn (matrices with 
to rows and n columns over S) componentwise. The monotonicity of addition 
and multiplication is preserved. Moreover, distributivity (|2"TT) also extends to 
matrices and vectors: 

A + sup.B'' = sup(A + B 1 - 1 ), A' • supB^ = sup(A' • B^), 
C + inf D» = inf(C + £>"), C ■ inf = inf(C" • D»). 

Here {B^}, {D^} are chains of matrices (ascending and descending, respec- 
tively), where {B^} is bounded from above. 

Indeed, the distributivity of addition is verified componentwise. Let us verify 
the sup-distributivity for multiplication. Let n be the number of columns of C . 
We have: 

n n 

(C • sup D»)ik = £ c' l3 ( sup d%) = sup ]T c'^ , (23) 
mga/" j=1 fj.aN' k . =1 

where A/" denotes a countable set and the last supremum is taken over all map- 
pings k from {1, . . . , n} to the natural numbers. The last equality is due to the 
scalar sup-distributivity. Now denote v := max™ =1 and observe that 

n n 

• £4^> (24) 

3=1 3 = 1 



since dj k are ascending chains. This implies that in the last supremum of (|23[) 
we can restrict maps n to identity, obtaining that 

n n 

sup^X-d^ = sup J2^ d % = sup(C'D) lk . (25) 

Thus the matrix sup-distributivity also holds. The inf-distributivity can be 
checked along the same lines replacing infimum by supremum, and ascending 
chains by descending chains. 

It can be checked that axioms A1-A3 and matrix distributivity (|2"21 provide 
sufficient ground for the proofs of Theorems 13.11 and 13.21 

For the alternative proof of Theorem 13.21 given in Remark 13.41 the system 
A1-A3 has to be modified. Note that the main part of the proof after (TT5)) does 
not need anything but the existence of sups of bounded ascending chains and the 
distributivity of addition over such sups. It is only the starting representation 
x = u + w, where u = A*b and w is the least vector w satisfying x = u + w, 
which may need more than that. 

We impose Al, A2 and the part of A3 asserting the distributivity of addition 
and multiplication with respect to sup's of ascending chains, which is needed 
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for Theorem 13. 11 that is, for u = A*b to be the least solution of x = Ax + b. We 
also require that there is at least one vector w satisfying x = u + w. This will 
hold if we impose 

A4. For all a, c € S such that a < c there is b £ S such that a + b = c. 

In addition, in order to get the least w satisfying x = u + w, we impose the 
distributivity of + with respect to arbitrary inf 's. Notice that in the case of an 
idempotent semiring we define the order canonically (a<b^a + b = b), and 
the axioms A2 and A4 are satisfied automatically. 

Now we consider some examples of semirings where Theorems 13.11 and 13.21 
are valid. In particular, axioms A1-A3 are satisfied for these examples. 

Examples 1,2: Max-times algebra, classical nonnegative algebra (see Pre- 
requisites). 

Example 3: Max-min algebra. Interval [0, 1] equipped with ab := 
min(a, b) and a + b := max(a, b). 

Example 4: Lukasiewicz algebra. Interval [0, 1] equipped ab := max(0, a+ 
6—1) and a + b := max(a, 6). 

Example 5: Distributive lattices. Recall that a lattice is a partially 
ordered set [5] where any two elements a, b have the least upper bound a V b := 
sup (a, 6) and the greatest lower bound a A b :— inf (a, b). A lattice is called 
distributive if the following laws hold: 

a V (6 Ac) = (a V 6) A (aVc), a A (b V c) = (a A 6) V (a A c) (26) 

When a lattice also has the lowest element e and the greatest element T, it can 
be turned into a semiring by setting ab := a A 6, a + b := a V 6, = e and 
1 = T. To ensure that the axioms Al and A3 hold, we require that the lattice 
is complete, i.e., that V Q a Q and f\pbp exist for all subsets {a a } and {bp} of the 
lattice, and that the distributivity can be extended: 

a V Apbp — Af}(a V bp), b A V a a a = V a (b A a a ). (27) 

Max-min algebra is a special case of this example. 

Example 6: Idempotent interval analysis. Suppose that a\,a,2 E S 
where S satisfies A1-A3, and consider the semiring of ordered pairs (01,02), 
where a\ < 02 and the operations of S are extended componentwise. This 
semiring, satisfying A1-A3, is the basis of idempotent interval analysis as intro- 
duced in [25] . 

Example 7: Extended order complete vector lattices. We can con- 
sider a semiring of all sequences (01, 02, • • ■) where € S and S satisfies A1-A3, 
with the operations extended componentwise. A generalization of Example 6 is 
then a semiring of all ordered sequences a\ < a<i < . . . where £ 5. 

Example 8: Semirings of functions. Further extension of Example 7 
to functions on a continuous domain is also evident (following [SJGS]). As an 
example of a subsemiring of functions satisfying A1-A3, one may consider convex 
functions on the real line, equipped with the operations of componentwise max 
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(as addition) and componentwise addition (as multiplication). In the spirit of 
max-plus semiring, we allow a function to take — oo values (which are absorbing). 
To verify A1-A3, recall that a function / is convex if the set {(x,t) | t > f(x)} 
is convex (providing connection to the well-known properties of convex sets) . In 
particular, the addition corresponds to the intersection of convex sets, and the 
multiplication corresponds to the Minkowski sum of convex sets. Note that the 
inf of descending chain of convex functions is also computed componentwise. 
As another example, we can consider a semiring, where an element is a class of 
functions on a continuous domain different only on a countable subset. Then, 
all countable sups or infs are well defined, since any two members of the class 
corresponding to such sup or inf will differ only on a countable subset, and 
axioms A1-A4 are verified componentwise, as above. 

Note that the Kleene star always converges in Examples 3-5. Moreover, it 
can be truncated for k > n, for anuxn matrix, so that A* = I + A+ . . . + A n ^ 1 , 
which follows from the optimal path interpretation, the finiteness of associated 
digraph, and because the matrix entries do not exceed 1. Hence A* b is well 
defined for any A, b, and x = Ax + b is always solvable, with the solution set 
described by Theorem 13.21 In Examples 6-8 the convergence of Kleene star 
should hold for the matrices corresponding to each component of the semiring 
(for the last "subexample" , excluding a countable subset of components). 

Finally we observe that theorems formally like 13.11 and 13.21 also hold in the 
case a linear operator leaving invariant a proper cone in R n , see [33], Theorem 
3.1, where an analog of Theorem 13.81 is also proved. 
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